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' Abstract. The spectral function of the scalar-isoscalar channel of the O(N) 

symmetric linear a model is studied in the broken symmetry phase. The inves- 
' tigation is based on the leading order evaluation of the self-energy in the limit 

. of large number of Goldstone bosons. We describe its temperature dependent 

variation in the whole low temperature phase. This variation closely reflects 
the trajectory of the scalar-isoscalar quasiparticle pole. In the model with no 
explicit chiral symmetry breaking we have studied near the critical point also 
Oh. the corresponding dynamical exponent. 
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1. Introduction 

The nature of the scalar meson sector of QCD, with particular emphasis on the 
broad a meson receives increasing attention [Si |^]. This sector describes the 
dynamics of the order parameter of chiral symmetry, which is the fundamental 
symmetry of quantum chromodynamics, realized in the broken symmetry phase. 
Therefore the response of this multiplet to the variation of the temperature or 
the baryonic density are directly related to the (partial) restoration of the chiral 
symmetry. 

We concentrate our study on the a-meson, therefore the 0(4) symmetric linear 
(T-model provides an adequate field theoretical framework. The finite tempera- 
ture/baryon density behaviour of the a meson has been studied extensively also 
with help of this model in the last 15 years [ |). Since in any realistic parametrisa- 
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tion the model is strongly (self )-coupled, there is no basis for the use of the simplest 
perturbative techniques. 

Improved (resummed) finite temperature perturbation theory was applied to 
the thermodynamics of the model [ and led to the conclusion that a becomes a 
narrow resonance with the increase of the temperature when its mass approaches the 
(temperature dependent) two-pion threshold. The narrowing is due to the smaller 
phase space available for the dominant a — > 2?t decay. This situation manifests 
itself in high energy heavy ion collisions by enhanced two-pion production near this 
threshold and the appearance of a high intensity narrow peak in the two-7 spectra 
superposed on the broad background arising from the ir° decay. 

One might remark, that results of the simplest perturbative treatments turn 
out to be sensitive to the details of the renormalisation prescription, do not reflect 
the so-called hybridisation phenomenon occuring between the a and some composite 
channels (see below), and finally, do not account correctly for the continuous nature 
of the chiral symmetry restoration. 

There was some notable progress in the non-perturbative lattice determination 
of the QCD spectral functions, which did not provide to date competitive results in 
the scalar-isoscalar channel [ §, q| ■ 

Our investigation relies on the large TV treatment of the thermodynamics of the 
linear sigma model, where TV refers to the number of the Goldstone bosons. The 
leading order approximation can be shown to avoid all problems of the usual pertur- 
bative treatment listed above. Additional credit for the use of large TV techniques is 
provided by its very successful application to dynamical critical phenomena occur- 
ing in superfluid helium [ |^, B |9| and by some very promising first attempts to go 
beyond Hartree-type approximations in exploring the mechanism of thermalisation 
for quantum fields [ . 

A more detailed report of our results obtained in the model with no explicit 0(4) 
symmetry breaking can be found in Ref. [ and a detailed investigation of the most 
realistically parametrised system with explicit symmetry breaking is presently near 
completion. Here we describe the more transparent (though phenomenologically 
less appealing) case when no explicit symmetry breaking is applied. 

2. Leading order large TV calculation of the propagators 

The Lagrangian of the TV-component scalar field theoretical model is parametrised 
for the large TV expansion in the following way: 

L = \\d^<f ~ mW] " ^(<T) V) 2 . (1) 

The direction of the symmetry breaking is chosen along the a = I direction and is 
expressed by the shift in the corresponding field component: 



(j> a -»• (VNQ + tp 1 ,^), i = 2,...,N. 



(2) 
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It is worth to note that in this paramctrisation 4> is related to the the phenomeno- 
logical f„ parameter at N = 4 by f n = 2$. 

The equation of state is determined by the vanishing of the coefficient of ip 1 in 
the shifted Lagrangian: 



m 



- i ^2 , ^ ll,„a\2\ 



+ —H^y\ . (3) 
6 ' 1 



To leading order in N the propagator of the transversal (Goldstone) modes receives 
the same tadpole contribution: 



A., _A_ 

6N 

which in view of (j^) leads to the masslessness of the Goldstone modes. 

The equation of state should be renormalised due to the ultraviolet divergence 
of (\/N)((ip a ) 2 }. To leading order in N one retains only the contribution from 
the massless Goldstone modes to the field fluctuations. The equation of state is 
rewritten with the renormalised couplings defined through the relations 

m 2 A 2 m\ 1 1 eA 2 1 

as 

< + ^$(T) 2 + ^§T 2 = 0. (6) 

This expression can be rewritten as an expression for the temperature variation of 
the order parameter relative to its T = value: 

® 2 (T) = $g - ^T 2 , (7) 

which yields for the critical temperature the prediction: T 2 = 12<& 2 = 3/ 2 « 
(161MeV) 2 . This is in very good agreement with the critical temperature deter- 
mined with lattice simulations [ 12 1. 

The propagator of the a field receives to this order additional contribution from 
the infinite chain of bubble diagrams: 

where b(p) stands for the analytic expression of a single bubble. After introducing 
the renormalised coupling, the renormalised, T — expression of the bubble is given 
by: 
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where Mq denotes the renormalisation scale. The T = cr-propagator looks like 

Xr _ n I 



(10) 



3 " 1 - Xr ln{-p 2 /M 2 )/96tt 2 ' 

In Fig.|l| we show the position of the T — complex tr-pole, determined from 
the equation G^ 1 = 0, as a function of Xr in units of f„. This pole is a physical 
resonance, since it has negative imaginary part. It is obvious from the figure that 
its real part cannot exceed ~ 3.5/ ff , while the imaginary part gradually increases 
with the growth of Xr. Choosing Xr — 310 we find M a jY a ~ 1, which is close to 
the phcnomenologically observed value, though the mass itself is too small. 

A pure imaginary pole (po = iM^) of the cr-propagator is found for every given 
T. The theory is meaningful only for momentum scales which lie considerably below 
the tachyonic frequency. This limitates the allowed Xr range from above. In Fig.|l| 
we also display the ratio Ml/Mu for several temperatures. The T-dependence is 
not very important therefore the limitation in the choice of Xr is nearly independent 
of T. 
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Fig. 1. The real and imaginary part of the T = pole mass of the a particle as a 
function of A_r. We also display the location of the pure imaginary tachyonic pole, 
M L . 



3. The scalar-isoscalar spectral function 



At finite temperature an additional term appears in the expression of b(p), which 
is defined in the upper halfplane of the complex po variable: 



1 



dx- 



1 



- 1 



1 



1 



Z + X 



Po_ 
2T' 



lmp >0. (11) 
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One can compute the real and imaginary parts of the complete expression of b(p) = 
bo(p) + bxip) when Impo — > +0. It determines the spectral function in the scalar- 
isoscalar channel: 

pOo,p) = — ImG a (p + iO,p). (12) 
Similarly to G a one can evaluate the propagator of the composite field (f> a (j) a — 

2. 

p 2 b(p)/6 + &(T)/3 
[P> p2(l~\ R b(p)/6)-\ R ^(T)/3- [ 1 

Obviously, its denominator is the same as that of G a , which reflects the phenomenon 
of hybridisation. 

In the two parts of Fig.|| we give the spectral functions of G a and of F for 
p = as functions of po for various temperatures. One observes that both T = 
curves are peaked at po values which agree with the real part of the cr-pole. The 
question we addressed was: can one interpret the shift of the maxima for increasing 
temperatures as a reflection of the finite temperature shift of the quasi-particlc pole? 

In order to answer this question one has to continue the propagators analyti- 
cally into the lower pq halfplane (onto the second Riemann sheet) and find the poles 
smoothly joining the location of the T = pole. The expression of the zero temper- 
ature piece b (p) is continued in a unique manner into this region. The harder task 
is the continuation of the finite temperature contribution. One can check that if one 
chooses br(po) to be analytic on the positive real axis, then the relation between 
b^ and its continuation into the lower halfplane b^ is the following: 

i i 

b^(po) =&£(po)- -. , /OTn r, lmp <0, Rep o >0. (14) 

47r exp(po/2J J — 1 

With help of this expression of b(po) valid on the second Riemann sheet one 
can find for every T the point where G r7 (po)~' 1 or _F(po) 1 vanishes. 

In Fig.|3| we display the imaginary part of the trajectory of the a-pole as a func- 
tion of the temperature. Its real part gradually decreases with increasing tempera- 
ture. The imaginary part slightly increases at the same time. At some T = Ti mag 
the pole position reaches the negative portion of the imaginary axis. At the same 
temperature a "mirror" pole touches the same point of the imaginary axis arriving 
from the left. The "collision" of the two poles leads to the birth of two purely 
imaginary poles, present for higher temperatures. One of them moves towards, the 
other one away from the origin along the imaginary axis. The narrowing of the a 
quasi-particle starts only for temperatures T > Ti mag and the same is true for the 
start of the "threshold enhancement" of the spectral function. 

The 0{N) symmetry is restored at the temperature for which the pole moving 
in this direction reaches the origin. In the temperature range when the poles are 
found on the negative imaginary axis, the spectral functions become cuspy. The 
closer moves the pole along the imaginary axis to the origin, the narrower the shape 
of the spectral function becomes. Finally for the temperature when one of the poles 
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Fig. 2. The spectral function of the elementary and of the quadratic composite 
scalar-isoscalar field as a function of the frequency po for different temperatures 
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Fig. 3. The imaginary part of the cr-polc in the 4 tn quadrant of the complex 
Po-plane including the negative imaginary axis. Depending on the details of the 
numerical algorithm used, the complex solution may continue into the up-going 
solution as well 



arrives to the origin one finds analytically that 

p( P(h T~T c )~S(p )/ Po , (15) 

where the <5-function singularity corresponds to the maximal threshold enhance- 
ment. 

From the point of view of the spectral function the whole temperature range 
Timag < T < T c can be called the threshold enhancement regime, since in this 
temperature interval the maximal signal comes from the direct neighbourhood of 
the imaginary axis. We believe that this is the generic behaviour, since there is 
no reason to expect that from every starting point in the po-plane the pole "flows" 
directly to the origin, without hitting the imaginary axis. 



4. Critical region 

In the vicinity of the critical point, where £ = T/87r<I> 2 (T) is the dominant length 
scale and the condition po, |p| << T is also fulfilled, one can derive an equation also 
for the soft modes with nonzero momentum: 
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Its solution in the approximation, when on the left hand side only the first term is 
retained exhibits the form of dynamical scaling [[?], H|: pa — |p| z /(|p|£) with z = 1. 
The second term provides the leading correction to scaling. In O(N) models, the 
dynamical exponent z = d/2 has been obtained for finite N on the basis of scaling 
and renormalisation group arguments, where in our case d = 3 [ [k| |lq] . For 
the correct interpretation of the situation it is important that there are two distinct 
hydrodynamical regions in the O(N) model for large N [|l3). In the true critical 
region z — d/2 is valid, in a precritical region z = 1 — 8Sd/Nd + 0(1/N 2 ), where 
Sd = 2/ir 2 for d = 3. The first region shrinks when N becomes large and completely 
disappears at N = oo. 



5. Conclusion 

The aim of this talk was to clarify the relationship of the temperature dependence 
of the spectral function characterizing the excitations of the chiral order parameter 
and the quasi-particle pole with the same quantum numbers (the cr-meson). The 
application of the leading order analysis in a large N approximation, where N is the 
number of Goldstone modes, was shown to avoid most of the conceptual problems 
related to the application of conventional perturbative approaches. The pole of the 
cr-propagator continued analytically into the lower po halfplane follows a trajectory, 
which describes satisfactorily the behaviour of the spectral function both near T = 
and in the critical region. The transition between the two regimes happens around 
Timag, the temperature when the pole becomes purely imaginary. 

The line of analysis presented here has been applied also to the model with 



explicit (9(4)-breaking external field [11, h6[, with parameters chosen the closest 



possible to the measured characteristics of the a — tt meson system. The behaviour 
found in this analysis proved generic also for the system with finite baryonic charge 
density. 

In view of the clear physical picture we found for the mechanism behind the 
threshold enhancement of the scalar-isoscalar spectral function, it is of interest to 
compute the next-to-leading corrections. This will enable us to see to what extent 
can one consider the present approach also quantitatively relevant. 
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